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Q^l Abstract. Classical novae are the results of surface thermonuclear explosions of 

hydrogen accreted by white dwarfs (WDs) from their low-mass main-sequence or red- 
giant binary companions. Chemical composition analysis of their ejecta shows that 
nova outbursts occur on both carbon-oxygen (CO) and more massive oxygen-neon 
^ , (ONe) WDs, and that there is cross-boundary mixing between the accreted envelope and 

underlying WD. We demonstrate that the state-of-the-art stellar evolution code MESA 
and post-processing nucleosynthesis tools of NuGrid can successfully be used for mod- 
eling of CO and ONe nova outbursts and nucleosynthesis. The convective boundary 
^ ' mixing (CBM) in our ID numerical simulations is implemented using a diffusion co- 

efficient that is exponentially decreasing with a distance below the bottom of the con- 
vective envelope. We show that this prescription produces maximum temperature evo- 
lution profiles and nucleosynthesis yields in good agreement with those obtained using 
the commonly adopted ID nova model in which the CBM is mimicked by assuming 
■^ . that the accreted envelope has been pre-mixed with WD's material. In a previous pa- 

^^ ' per, we have found that ^He can be produced in situ in solar-composition envelopes 

^ ; accreted with slow rates (M < 10-^°Mo/yr) by cold (T^wd < 10^ K) CO WDs, and 

that convection is triggered by ^He burning before the nova outburst in this case. Here, 

we confirm this result for ONe novae. Additionally, we find that the interplay between 

, , , the ^He production and destruction in the solar-composition envelope accreted with an 

r> ■ intermediate rate, e.g. M - lO^'^Mo/yr, by the I.ISMq ONe WD with a relatively 

3 ' high initial central temperature, e.g. Two = 15 x 10^ K, leads to the formation of a thick 

- - - radiative buff'er zone that separates the bottom of the convective envelope from the WD 

surface. 



o 



1. Introduction 

MESA is a collection of Fortran-95 Modules for Experiments in Stellar Astrophysics. 
Its main module star can be used for ID st ellar evolution s imulations of almost any 
kind (numerous examples are provided by Pa xton et al.l(l201ll l2013) and on the MESA 
website http : //m esa . sourcef orge . n et). For instance, star can compute with- 
out any interruption the evolution of a solar-type star from the pre-MS phase through 
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the He-core flash and thermal pulses on the asymptotic giant branch towards the WD 
cooling. Only a few codes can do such computations. Other MESA modules provide 
star with state-of-the-art numerical algorithms, e.g. for adaptive mesh refinement and 
timestep control, atmospheric boundary conditions, and modem input physics (opaci- 
ties, equation of state, nuclear reaction rates, etc.) 

The nucleosynthesis calculations presented in this work are made by the post- 



proce ssing code PPN (developed inside the NuGrid research framework, iHerwig et al.l 



20081) . where the input data of stellar structure (e.g., from MESA models) are processed 
using a comprehensive list of nuclear species and reaction rates. The network can in- 
clude more than 5000 species, between H and Bi, and more than 50000 nuclear reac- 
tions. A self-controlled dynamical network defines the actual numbers of species and 
reactions considered in calculations, based on the strength of nucleosynthesis fluxes. 
For the complete post-processing of the full MESA nova models, we use the multi-zone 
parallel frame of PPN (MPPNP), whereas for the simulations based on a single-zone 
trajectory we use the serial frame (SPPN). The two frames use the same nuclear physics 
library and the same package to solve the nucleosynthesis equations (for details, go to 
http : //www . nugridstars . org). 

We have prepared shell scripts, template MESA inlist files, and a large number 
of initial WD models aiming to combine MESA and NuGrid into an easy-to-use Nova 
Framework. This new research tool can model nova outbursts and nucleosynthesis oc- 
curring on CO and ONe WDs using up-to-date input physics and a specified nuclear 
network for a relatively large number (up to a few thousands) of mass zones cover- 
ing both the WD and its accreted envelope. Nova simulation results can be analyzed 
using animations and a variety of plots produced with NuGrid 's Python visualization 
scripts. This article reports the new results obtained with the Nova Framework and 
demonstrates some of its capabilities. 

2. One-Dimensional Nova Models With Convective Boundary Mixing (CBM) 

To be consistent with observations, a realistic nova model should assume that the ac- 
creted material has solar composition (for the solar metallicity) and that it is mixed with 
WD's material before or during its thermonuclear runaway. Recent two- and three- 
dimensional nuclear-hydrodynamic simulations of a nova outburst have shown that a 
possible mechanism of this mixing are the hydrodynamic instabilities and shear-flow 
turbulence induced by steep horizont al velocity gradients at the bo ttom of the convec- 
tion zone triggered by the runaway (ICasanova et al.ll2010l l2011al Jbh. These hydrody- 



namic processes associated with the convective boundary lead to convective boundary 
mixing (CBM) at the base of the accreted envelope into the outer layers of the WD. 
As a result, CO-rich (or ONe-rich) material is dredged-up during the runaway. In our 
ID MESA simulations of nova outbursts, we use a simple CBM model that treats the 
time-dependent mixing as a diffusion process. The model approximates the rate of mix- 
ing by an exponentially decreasing function of a distance from the formal convective 
boundary, 

OcBM^A)expl -— — , (1) 

where Hp is the pressure scale height, and Do is a diffusion coefficient, calculated using 
a mixing-length theory, that describes convective mixing at the radius ro close to the 
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boundary. In this model / is a free parameter, for which we use t he same value f = 
fnnva - 0.004 that we used in our CO nova simulations in Paper I ("Denissenkov et al" 



l 2013h . For a physical motivation to model the CBM with the prescription ([T]), see 
Paper I. Our MESA nova models with the CBM have been post-processed with the 
MPPNP code, and the resulting final abundances have been compared with those ob- 
tained for the corresponding nova models without CBM but with the 50% pre-mixed 
accreted envelopes. The comparison between the final abundances shows a very good 
agreement for both the CO and ONe nova models (the left panel of Fig.[T]and Fig. 1 in 
IPenissenkov et al.ll201Z : we disavow our conclusion made in that paper that it was not 
true for CO novae). This gives a support to the widely-adopted ID nova model in which 
the CBM is mimicked by assuming that the accreted envelope has been pre-mixed with 
WD's material. 



Mass range [1.150005-1.150025], Cycle 1200 



g. 10 

o lO' 
S lO' 

01 

^ ion 

-a 

111 
< 

.lio 



: 10- 



tii 



•j^; 



i 



15 20 25 30 
Mass Number 




Figure 1. Left panel: A comparison of final solar-scaled abundances of stable 
isotopes from two of our 1.15 M© CO nova simulations with Two = 12 MK and M = 
2 X lO^'^Mo/yr. The red circles are the results obtained with the 50% pre-mixed 
initial abundances in the accreted envelope, while the blue squares represent the 
case with the solar-composition accreted material and convective boundary mixing 
(CBM) modeled using the same method and CBM parameter /nova - 0.004 that 
we used in the 1.2 Mq CO nova model in Paper I. Right panel: Tniax-trajectories 
from three of our 1 .3 Mq 50% pre-mixed ONe nova models computed for the same 
accretion rate, M - 2x 10"'" Mg/yr, but different WD's initial centtal temperatures. 
The i nitial isotope abundances are those of the Barcelona group jJose et alj|1999l 
l200Th . 



3. Comparison of the Results of SPPN and MPPNP Simulations 



Single-zone post-processing calculations of nova nucleosynthesis with the NuGrid code 
SPPN are done for a mass zone in a nova envelope where the temperature profile has its 
maximum, T = r^ax (usually, near the core-envelope interface). For this, the tempo- 
ral variations of Tmax and its corresponding density, called "trajectories", are extracted 
from MESA models (e.g., the right panel of Fig. [T}. The SPPN code runs much faster 
than MPPNP, therefore it can be used for a comprehensive numerical analysis of param- 
eter space, when the relevant isotope abundances and reaction rates are varied within 
their observationally and experimentally constrained limits. A comparison between 
the final abundance distributions in Fig. |2] shows a qualitative agreement between the 
results from the single-zone rmax-trajectories (the left panel) and from their correspond- 
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Figure 2. Left panel: the final solar-scaled abundances of stable isotopes cal- 
culated with the SPPN code using the trajectories from the right panel of Fig. [T] 
that correspond to the three different WD's initial central temperatures: Two - 12 
MK (diamonds), 15 MK (squares), and 20 MK (circles). Right panel: the same fi- 
nal abundances, but calculated with the MPPNP code using the corresponding full 
MESA models and averaged over the indicated mass range in the expanding nova 
envelope. 



ing multi-zone complete models (the right panel). Indeed, the abundance patterns for 
different groups of isotopes look similar. The main quantitative difference is the more 
efficient production toward heavier species in the first case, where O is more depleted, 
while Ar and K isotope^ are more efficiently made. This difference is caused by the fact 
that in the multi-zone post-processing simulations convective mixing reduces an "aver- 
age temperature" at which the nucleosynthesis occurs and it also constantly replenishes 
the hydrogen fuel burnt at the base of the envelope by bringing it from its outer parts, 
where T <c Tmax- In general, the results in Fig. |2] show how the increase of the r^ax 
peak value (the right panel of Fig. \Q changes the relative abundance distribution. We 
conclude that SPPN nucleosynthesis simulations with nova Tmax-trajectories provide a 
proper qualitative indication of the behavior of their corresponding complete MPPNP 
simulations. However, they can be used only as a diagnostic for nova nucleosynthesis, 
e.g. in reaction rate sensitivity studies. 



4. Effects Caused by ^He Burning 

In Paper I, we have found, for the first time, that in the extreme case of a very cold CO 
WD, e.g. with T\f^D = 7 MK, accreting solar-composition material with a very low rate, 
say M = lO~^'M0/yr, the incomplete pp I chain reactions lead to the in situ synthesis of 
^He in a slope adjacent to the base of the accreted envelope and that the ignition of this 
^He triggers convection before the major nova outburst. The ^He burning continues at 
a relatively low temperature, T ~ 30 MK, approximately until its abundance is reduced 
below the solar value, only after that the major nova outburst ensues triggered by the 
reaction ^^C(p,y)^^N. Although this is an interesting variation of the nova scenario^] 
because the CO enrichment of the accreted envelope is produced by the ^He-driven 



'The Barcelona composition specifies the initial abundances only for isotopes lighter than Ca. 
-Initially, ^^He was considered as the most likely isotope to trigger a nova outburst bv lSchatzmanI jl951h . 
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CBM in this case, its observational frequency is expected to be very low (Paper I) . 
Later, we have confiiTned this result for ONe novae (Fig. 2 in iDenissenkov et alj2012h . 



Furthermore, we have found that in the 1.15 M© ONe WD with a relatively high initial 
central temperature, Twd - 15 MK, accreting the solar-composition material with the 
intermediate rate, M = 10^'" Mo/yr, the interplay between the ^He production and 
destruction in the vicinity of the base of the accreted envelope first leads to a shift 
of Tmax away from the core-envelope interface followed by the formation of a thick 
radiative buffer zone th at separates the bottom of the convective envelope from the 



WD surface (Fig. 3 in IDenissenkov et al.ll2012l) . This result is obtained only for the 



case when an ONe WD accretes solar-composition material and it almost disappears 
in the models with the pre-mixed accreted envelopes. Given that the formation of the 
radiative buffer zone is revealed in the computations with the more realistic nova model 
parameters, such cases can probably be observed. We defer a more detailed study of 
this peculiar case and its possible consequences to a future work. In particular, it would 
be interesting to see if the CBM can cross the buffer zone and reach the WD. 

5. Conclusion 

We have created the Nova Framework that allows to simulate the accretion of H-rich 
material onto a white dwarf (WD) leading, for a suitable set of initial parameters, to a 
nova outburst, and to post-process its accompanying nucleosynthesis. The Nova Frame- 
work combines the state-of-the-art stellar evolution code MESA and post-processing 
nucleosynthesis tools of NuGrid. It includes a number of CO and ONe WD models 
with different masses and central temperatures (luminosities) that can be used in simu- 
lations of nova outbursts. The use of the Nova Framework is facilitated by a number of 
shell scripts that do the routine job necessary to coordinate the operation of the MESA 
and NuGrid codes. 

Within the Nova Framework project, we have provided a large set of nucleosyn- 
thesis calculations for nearly 50 CO and ONe nova models at the solar metallicity 
(Z = 0.02) so far. To verify our calculations, we have compared their results with 
those published in the literature for similar nova models. The comparison shows a very 
good qualitative agreement. 
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